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A general expression is discussed for pion probability distributions coming from relativistic heavy 
ion collisions. The general expression contains as limits: 1) The disoriented chiral condensate (DCC), 
2) the negative binomial distribution and Pearson type III distribution, 3) a binomial or Gaussian 
result, 4) and a Poisson distribution. This general expression approximates other distributions such 
as a signal to noise laser distribution. Similarities and differences of the DCC distribution with these 
other distribution are studied. A connection with the theory of disordered systems will be discussed 
which include spin-glasses, randomly broken objects, random and chaotic maps. 
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The purpose of this paper is to discuss a general expression for the pion probability distribution which may be used 

to analyze pions coming from relativistic heavy ion collisions. The study of pions is of current interest for several 

reasons. First pions are the main component of the produced particles coming from such collisions. Several thousand 

>" , pions are now observed at CERN SPS experiments and this number may go up by a factor of 10 at RHIC energies. 

[--[- ■ Secondly, the behavior of pions may signal the formation of the quark-gluon plasma as, for example, in the disoriented 

' chiral condensate (DCC) picture K-M ■ Thirdly, the fluctuations in pions have been discussed in terms of intermittency 

O^ • and fractal structure [[7| coming from non-Poissonian effects. A distribution which has been used to discuss these 

, I phenomena is the negative binomial (NB) distribution mM- The NB distribution also has an important feature 

►^ , known as Koba-Niclscn-Olessen (KNO) scaling [|lO[ which the Poisson distribution lacks. Fourthly, the Bose-Einstein 

^^ correlations amongst pions is an important property used in Hanbury-Brown-Twiss (HBT) experiments, and such 

p\| , correlations have also been proposed for the formation of a pion laser ||ll|,|l^ . Bose-Einstcin condensation of atoms 

^^ ' in a laser trap is a very recently observed phenomena in another area of physics. Distributions which extrapolate 

C^ ' between Poisson and Bose-Einstein and negative binomial have been developed, such as the signal to noise model of 

^^ Glauber-Lach. This S/N model was originally developed in quantum optics, but has also used for particle production 

and a review can be found in [ [l3[ . The expression to be developed extends the range of extrapolation by including 

not only these distribution but others. Specifically, a general expression will be developed which contains as special 

i-^ limiting cases the DCC distribution, the NB distribution, a binomial distribution p4| or Gaussian like distribution, 

I , and the Poisson distribution. The general distribution also approximates the S/N model. A connection with the 

fi ' theory of disordered systems p5|Jlfl] w ill also be discussed further. These disordered systems included spin-glasses 

^ |0,|8[, randomly broken objects [[IQpcl ], random permutations [2^, random maps |^,|9[, and chaotic maps p2[ . 

^ For example, the disoriented chiral condensate gives a probability distribution which comes from a random direction 

►^ [ of the isospin vector and a connection of the DCC state with the theory of disordered systems will be noted. The 

• ^ . probability distribution is simply given by 

X- 

■^-' P,iN,.,j)^^^j g-^-^— ^ (1) 

where fc is a random variable which can take on values fc = 0, 1, 2, ..., N, and x, 7 are parameters. The Beta 
function that appears in Eq.(n|) is Beta[ci;,2;] = r(uj)T(z)/r(uj + z) where r(z) is a gamma function. In probability 
theory, Pk{N,x,j) is known as a Polya distribution ||2J]. This distribution was used in Ref. (25| to describe cluster 
yields coming from high energy heavy ion collisions and a connection of it with the theory of disordered systems was 
discussed briefly in Ref. |l^ . Its usefulness as a model for pion probability distribution will be developed. 

Some interesting hmits of Eq.dll) are as follows. Setting x = 1, j — 1/2 in Eq.P), the Pk{N, x, 7) = Pk{N, 1, 1/2) 
is given by 

(7V!)222^(2/c)! 1 

P,(A^, 1,1/2)--^-^—-^-^^-^=. (2) 



{2N + l)l{k\2''y 

Eq.(0) appears in pionic yields from the disoriented chiral condensate model of a QCD phase transition. Specifically 
2k = riQ is the number of neutral pions and 2A^ — uq + n+ + n_ is the total number of pions with n+ = n_ being 
the number of positive or negative pions. In obtaining l/2\/Nk, Stirling's approximation was used. The probability 
distribution for the neutral to total pion yield, R^ — no/(no + ^1+ + 'i-)j is then PiRs) ~ 1/2a/R^. A general 
expression for the fluctuation of Eq. (nh is |l5[ 



<'^>-<'>'-JT7TT«0 + 



f) (3) 



where (fc) = (7/(2; + 7)) iV. For a; = 1, 7 = 1/2, (fc> = (1/3)A^ or (no) = 2Ar/3 = (n+) = (n_). Thus (ng) - {naf = 
(4/5) (no) (1 + (no)) and (nl) - {n+f = (1/5) (n+) (1 + (n+)) = (l/4)((n2) - (no)^). The probabUity distribution of 
j TT+'s or 7r_'s is 

C^h222j (2(N - j)V 1 

^.l^v,i,i/z; ^2N + l)\{{N-3W 2^N{N-j) ^' 



The Pj{N, 1, 1/2) increases with j and behaves asymptoticaUy hke the arcsine distribution 1/ti^ j{N — j) for j ^ N 
when A^ is very large. This behavior is totaUy different than the ttq behavior simply obtained from Eq.(g). As wiU be 
shown the tto behavior can be approximated by other distributions, but not this 7r_|_ or 7r_ behavior which is coupled 
to the ttq distribution. The arcsine distribution is also a special case of a distribution related to Eq.(n]) to be discussed 
below. The constraint 2N = no + n+ +n_ represents a major difference between the DCC model and a NB description 
and other probability distributions usually applied to pionic yields. 

A binomial limit of Eq.(|l|) is obtained in the limit x -^ 00, 7 — *■ 00. If a; = 27 so that (no) = {n+) = (n_) = ■i(2iV), 

then Pk{N) = ( ^ ) {^t (1)^ ''■ ^1^° ('^^) - ^^f = (2/3) {k) from Eq.(|I|) and since no = 2fc: (ng) - (no)^ == 

(4/3) (no). For large iV, the binomial distribution is approximately a Gaussian. The DCC distribution for a large 
number of domains is also a binomial. In the limit in which each domain has fc = or 1 or no = 0, 2 tto's so that 
the number of domains equals N , then the Pk{N) is also given by the binomial result quoted above. Moreover, the 
DCC distribution has the following central limit theorem feature. Let m equal the number of domains and let A'^i 
be the maximum fci in any one domain, i.e., fci — 0, 1, 2, ..., iVi, with noi = 0, 2, ..., 2Ni the random number of 
TTo's coming from each domain. The fluctuation in one domain is (fci) — (ki) = af and is given by Eq.(Bf) with 
x = 1, 7 = 1/2, or using (ni) = iVi/3: (kf) - {kif = (2/15)7Vi(l + 27Vi/3). For m DCC cells, the total variance is 
al^ = mal = (2/15)7V(l + 27V/3m). When m = iV, cr^ = {2/9)N which is the binomial result. 

A negative binomial limit can also be obtained from Eq.dl^) in the limit N^oox-^(X),p — x/N and p — p/{l + p). 
Then 



mp,7)=['^t-' 



)p^(l-p)' 
7 + k-l\ 1 / {k)h 



k ) {i + {k)hy\i + {k)h 



(5) 



where {k) = (7/(2; + 7)) iV = 7/p and (fc^) - {kf = {k) (1 + {k) h) from Eq.(|). The Pk of Eq.(|) is the NB 
distribution with 7 now the NB parameter. The NB distribution is of current interest since this distribution fits many 
experimental observations as, for example, in the multiplicity distribution of hadrons produced in e^e~ collisions at 
LEP and PEP-PETRA and in hadron and heavy ion collisions [26|. Van Hove and Giovannini [^,^ have discussed a 
clan model which gives a NB distribution. Carruthers and Shih |12] summarized various mechanisms that lead to NB 
distributions and show that this distribution can be put into a correspondence with self-similar Cantor sets or fractal 
structures. The NB distribution has been used to discuss issues related to intermittency [0j|]. Taking 7 = 1/2, the 
same value of 7 as in the DCC distribution, then 

Pji/,)- (2fc)! ( 1 V^7 nk) \ 

The factor (2A;)!/2'^(A:!)^ also appears in the DCC distribution and its Stirling limit is l/y/nk which gives the 

fc 



1/y/R^ behavior associated with the DCC state. The factor (2 (fc) /(I + 2 (fc))) gives a geometric series decrease 
or exponential decrease in Pk- Writing (2 (fc) /(I + 2 (fc))) = e^°'^ with a — ln(l + 1/2 (fc)) « 1/2 (fc) results in 

Pfe(l/2) « {l/^2Tik{k)\ e-'^/^W. For general 7, 

Pkil) « (fc^-Vr(7)) (7/ {k)y e-^'^/W (7) 

which is a Peason type III distribution. 



The Poisson limit of the NB is reahzed when 7 ^ (fc). Both the NB and DCC have larger than Poisson fluctuations. 
Fluctuations of pions larger than Poisson are to be expected just on the basis of Bose-Einstein statistics, with the 
Poisson limit obtained from Maxwell-Boltzmann statistics. A 20% enhancement in fluctuations above those of Poisson 
statistics occurs in both hydrodynamic and thermodynamic models from Bose-Einstein correlations |27| ] . An important 
feature of the NB and Peason type distribution is that (fc) Pk{j) is a universal function of the scaled variable k/ (k) 
in the limit k ^ 00, (k) ^ 00, which is known as KNO scaling ijl^l. For example the Peason Type III is simply 

1 /2 

kPkil) = (7^/r(7))(fc/(fc))'^"^e-'^'''/<'=> and the DCC distribution is (k) Pk = (1/2V3) (^] . However, the 
Poisson distribution does not obey KNO scaling. 

Other distribution have been developed which extrapolate between Poisson and NB and have been used to describe 
the multiplicity distribution of produced particles. These include the Glauber-Lach signal S to noise NL model 
| fi3[ with S the coherent signal level (Poisson emitter) and NL the thermal Bose-Einstein noise level. In Biyajima's 
generalization of this model [[l3| , P8[ , the probability of k particles is given by a distribution which is called a laser 
distribution: 



Pik) 



[NL/aY 



[l + NL/a) 



k-\-cx. 



^-S/(l+i^)^a-l 



"- NL 



1 + ^ 



(8) 



with a = 1 the Glauber-Lach model and i^^^ is a generalized Laguerre function. The (fc) = NL + S and (fc-^) — (fc) — 
{k)+NL{2S+NL)/a. The Poisson limit is achieved when NL -^ 0, while a NB limit follows for S ^ 0. By comparison, 
the distribution of Eq.(|^) also contains these two limits as well as the DCC distribution and other distributions already 
mentioned. Fig.l illustrates the similarity of Eq.dlf) for the DCC choice of a; = 1, 7 = 1/2 and for N = 300 with the 
distribution of Eq.(|) with a = 1/2, S = 62, NL = 38. Both distributions have (fc) = 100. In the DCC choice k is 
restricted to < fc < iV while for signal to noise models < fc < 00. The figure shows the interval < fc < 300 only. 
As noted, a major difference between the DCC results and these other models is the constraint 2A^ = ttq -I- 7r+ -f 7r_ . 
This leads to totally different behaviors in the ttq channel and the 7r_|_, 7r_ channel. This result suggests that event- 
by-event data should be investigated not only for the fall off in ttq to total yield as 1/2-^^ but also the rise in the 
7r+ or 7r_ yields as given by Eq.ffl). Also shown in this figure is a NB distribution with 7 = 1/2 and (fc) = 100. If the 
fluctuations of the two distributions, Eq.([I]) and Eq.(||) are equated then a, 7 are connected by 7 = a/(l — b^) where 
b = S/{NL + S) is the fractional signal level. Very large Bose-Einstein enhancements to Poisson results can appear 
in other models such as the pion laser model of Ref. 
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FIG. 1. The distribution of Eq.(h|) for the DCC choice a:: = 1, 7 = 1/2 (solid line) compared with the negative binomial 
distribution with 7 = 1/2 and (k) — 100 (dashed line) and a signal to noise model with 5* = 62, NL — 38 and a — 1/2 (dotted 
line). The three distributions have similar features. 



When Eq.(^ is rewritten as 



Pk{N,x,-f) 



N 



w {1 — w) u{w,x,j)dw 



(9) 



important functions as M(w,a;, 7) = (Beta(x, 7)) w^ ^{l — w)^ ^ and f{w,x,j) — w ^u{w) emerge. The u{w,x,j) 
is also the limit of Pk{N, x, 7) as iV ^ 00, fc — > 00, k/N —^ w. The Stirling limit of Eq.(n]) corresponds to u{'w, 1, 1/2). 



The distribution of Eq.(S) and u{'w, x, 7) and its associated /(w, x, 7) play a prominant role in the theory of disordered 
systems such as spin-glasses, randomly broken objects, random and chaotic maps, etc. Table 1 summarise some of 
these connections. The choice of x and 7 varies for the different areas as listed in the table. The last column gives the 
quantity described by Eq.(|9|). In spin-glass models, random hamiltonians based on an Ising interaction Ja^ • aj (where 
J = ±|J| is chosen randomly) are used to calculate rugged free-energy landscapes. The function f(w,x,"f) and its 
related probability function u{'w,x,"f) give the distribution of well depths in the free energy landscape. For clusters 
and breaking processes, with initial size A into fragments of size k — 1, 2, ..., A with Uk ~ the number of fragments of 
size fc, f{w,x,"f) is related to Uk by Auk — > f{w,x,j) in the limit fc, A — > 00 with w — k/A. At a; = 1, n^ '--^ 1/fc^, a 
power law distribution with Fisher exponent r = 2 — 7. In the case of random permutations in Table 1, the n^ of the 
cluster example just given now becomes the number of cycles of length k, when the symmetric group of permutations 
is represented by its cycle class structure. A cycle class description for pionic distributions was also described in Rcf. 
p8[ . The random map example of Table 1 is also a prototype model for disordered systems [Q. Here n^ is the 
number of attractors with k elements. These attractors are limit cycles produced in a random mapping of A points 
into itself (allowing multiple points to go into one point). The f{w, x, 7) and u{w, x, 7) = wf{w, x, 7) follow as in the 
cluster case above. Chaotic maps, such as the tent and quadratic map |£3|,^, have invariant distributions which are 
given by arcsine laws. The u{w, x, 7) for x — j = 1/2 is the arcsine distribution u{w, x, 7) = n/ (w(l — wj) 

TABLE I. Various disordered systems and values for x, 7 described by Eq.(pl). 



Area 


X 


7 


Quantity Described 


Spin Glass 0,|l|| 




1-1-7=1 


Well depths from random hamiltonians 


Disoriented Chiral Condensate - DCC [ 1 -|5 


1 


1/2 


Pion probability distribution 


Randomly Broken Objects [|6|jl|,|^ 


X 


1 


Fragment sizes 


Random Permutations [gl| 


1 


1 


Cycle length 


Random Maps jlf ,1§| 


1/2 


1 


Limit cycle distribution 


Chaotic Maps; Tent and Quadratic |2^, ^ 


1/2 


1/2 


Invariant distribution 


Power Laws with Fisher Exponent | L5| 

r < 2 and r = 2 - 7 


1 


7 


Cluster or droplet distributions 



To summarize and conclude, this paper contained a general expression for pions coming from relativistic heavy 
ion collisions. Various limits of this general distribution contain frequently used pion distribution. One important 
example is the disordered chiral condensate distribution which is of current interest because of its use as a possible 
signal of the quark-gluon phase transition. Another important limit gives rise to the negative binomial distribution 
which appears in discussions of intermittency, KNO scaling, non-Poissonian fluctuations. Other limits include a 
binomial limit, or Gaussian like limit, and a Poisson limit. Distributions which extrapolate between Poisson and 
negative binomial such as the Glauber-Lach and laser distribution have been developed previously and this paper 
extends the range of limiting cases to include many other distributions. Comparisons are made between the DCC, 
NB distribution, and laser distribution, with each having large fluctuations. Similarities and differences which may 
be important in understanding the quark-gluon phase transition are discussed. For example, the DCC distribution 
has a probability distribution of ttq's which falls a l/vA:, Eq.(0), coupled with a 7r_|_ or 7r_ distribution which rises as 
1/y/N — fc, Eq.(|j), a,s j ^ N for j tt+'s or 7r_'s. The ttq distribution can be approximated by a standard negative 
binomial model or a signal to noise laser model as shown in Fig.l. Thus, this unique coupling of the ttq and 7r+ or 
7r_ channel characterizes the DCC state more than just the often quoted l/y/Rz or 1/vk behavior in the probability 
distribution of ttq's. Heavy ion pion production data may also be analyzed in terms of the general distribution. The 
connection of this general distribution with the theory of disordered systems is discussed. These connections included 
spin glasses, randomly broken objects, random permutations, random maps and chaotic maps. 
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